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r— I ■ Abstract 

^ i' We analyze general enough models of repeated indirect measure- 

ments in which a quantum system interacts repeatedly with randomly 
chosen probes on which Von Neumann direct measurements are per- 
formed. We prove, under suitable hypotheses, that the system state 
probability distribution converges after a large number of repeated 
indirect measurements, in a way compatible with quantum wave func- 
tion collapse. Similarly a modified version of the system density matrix 
^ , converges. We show that the convergence is exponential with a rate 

' given by some relevant mean relative entropies. We also prove that, 

, under appropriate rescaling of the system and probe interactions, the 

state probability distribution and the system density matrix are so- 
lutions of stochastic differential equations modeling continuous-time 
quantum measurements. We analyze the large time convergence of 
these continuous-time processes and prove convergence. 
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1 Introduction 



^ ■ Repeated indirect quantum measurements aim at getting (partial) informa- 

■ ■' tion on quantum systems with minimal impact on it. A possibihty consists 

in repeating non-demolition measurements (QND). At each step, one lets the 
quantum system under study interact with another quantum system, called 
the probe, and performs a Von Neumann measurement [1] on this probe. In- 
formation on the quantum system is gained through intrication between the 
probe and the quantum system. If one is aiming at progressively measuring 
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a quantum observable, one has to make sure that a system prepared in one 
eigenstate of this observable remains in it after a cycle of intrication and 
direct measurement on the probe, and that the set of stable states, called 
pointer states, forms an orthonormal basis of the system Hilbert space. The 
experiment of ref.[2], in which repeated QND measurements is used to fix 
and measure the number of photons in a cavity without destroying them, 
illustrates this strategy. 

Repeated indirect measurements were studied in ref.[T]. There, the dis- 
cussion was limited to QND measurements consisting of identical probes, in- 
teractions and measurements on the probe, and assuming a non-degeneracy 
condition. In the present article, we extend these results to cases where dif- 
ferent indirect measurements (probes, interactions and direct measurements 
on probes) are used. We also study the degenerate case. 

Let Qn{') be the pointer state distribution after the n*^ indirect measure- 
ment, that is Qn{ct) is the probability to find the system in the state a after 
n steps (a labels the pointer states), Qnict) = 1- As explain in section [3l 
each cycle of indirect measurement updates the distribution through Bayes' 
law. The analysis of the distribution is reformulated as a problem in classical 
probability theory (with no quantum interference). 

We shall prove that this sequence of distributions converges at large n, 
that is after a large - strictly speaking, infinite - number of QND measure- 
ments. If a non-degeneracy assumption is verified, the limit distribution is 
Q{a) = 6r,a for some random limit pointer state T. This reflects the collapse 
of the system wave function. The convergence is exponential, 

(5n(a) ~ const. exp(— n5(T|a)), for large n, a 7^ T, 

with rate given by an appropriate relative entropy S{T\a) deflned in eq.Q. 
In probabilistic terms, the limit Q{a) possesses a natural interpretation as a 
Radon-Nikodym derivative and, Qn{oi) = '^{Q{a)\J-n) is a closed martingale 
with respect to an appropriate flltration, see section [5j As a consequence, 
we show that the expectations conditioned on the limit pointer state are 
identical to expectations starting from this same pointer state. That is: 
Ex( • ) = IE( • 1^) where A is the tail c-algebra, the smallest one making the 
limit distribution measurable. See below for a precise definition of Ex( • )■ 

Convergence of Qn{ot) is also studied when the non-degeneracy hypothe- 
sis is not fulfilled. In this degenerate case, the limit pointer state distribution 
vanishes outside a random finite set of pointer states. The quantum system 
density matrix, when properly modified, also converges in the limit of infinite 
number of QND measurements. The limit density matrix then coincides with 
that predicted for degenerate Von Neumann measurements Uj, see section [5j 

Of course repeated indirect measurements have already been studied in 
the physics literature, mostly through time continuous measurement for- 
malisms - as far as we know, little was done on the discrete setting as we 
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do in the present paper. E. B. Davies[6] probably made the first rigorous 
approach to time continuous quantum measurement. This was later studied 
by N. Gisin[7] and L. Diosi[8j using the non linear Schrodinger equation. 
Simultaneously, A. Barchielli and V. P. Belavkin derived the equations gov- 
erning continuous measurements in terms of instruments [9]. They derived 
jump equations which, when properly rescaled, are equivalent to diffusive 
equations for continuous measurements. Another approach uses quantum 
stochastic differential equations and quantum filtering theory to obtain the 
so-called Belavkin equations \10\ \TT\ \T7[ [T3| [T^ and |15) . More recently, 
C. Pellegrini derived Belavkin equations for continuous time measurements 
|20| [2H [22] using discrete repeated indirect measurement models. The prob- 
lem of convergence of quantum density matrix has also been analyzed within 
the time continuous measurement framework. In refs.|161[T7]. V. P. Belavkin 
showed the convergence of mixed states toward pure states. A derivation of 
wave function collapse from the non-linear stochastic Schrodinger equation 
has been presented in refs.|18| 119) . It makes use of martingale theory as we 
do in the present paper. 

In the following we also connect our discrete model to the time con- 
tinuous measurement formalism. Taking the time continuous limit requires 
rescaling appropriately the interaction between the quantum system under 
study and the probes. In that sense, the time continuous model we consider 
is close to that of ref. |21| but our proofs are different and slightly more gen- 
eral. Our derivation is based on the convergence of some discrete counting 
processes - related to the number of occurrences of outputs in the successive 
indirect measurements - toward a time continuous Gaussian process. Under 
appropriate hypotheses, spelled out in section [6l the pointer state distri- 
bution satisfies a random diffusive stochastic equation driven by Gaussian 
processes. Suppose that at each step the probe system is randomly selected 
(independently of the past history and with time independent probability, 
for simplicity) among a finite set O whose elements are indexed by o G O 
and that the output measurements on the probe can take finite number of 
values indexed by i G spec(o). Then, the pointer state probabilities Qt{a) 
are time continuous martingales (with respect to an appropriate filtration) 
whose evolutions are governed by the non linear stochastic equations: 

dQtia) = gi(a)^ (r(")(i|a) - {T^°\i})t)dXtio,i) 

(o,i) 

where Xt{o,i) are some centered Gaussian processes, T^°\i\a) are coding 
for the probability of output probe measurement i within the probe sys- 
tem o conditional on the quantum system be prepared in the state a and 
{T^°\i))t = ^^T^"\i\a)Qt{ct)- The pointer state distribution again con- 
verges as a finite-dimensionnal bounded vector martingale. Under non- 
degeneracy assumptions, the limit distribution is again Q{a) = 6r,a and 
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the convergence is still exponential with a rate given by the scaling limit of 
the mean relative entropy. 

These results extend to the system density matrix. In the time continuous 
scaling limit, the system density matrix is a solution of a diffusive Belavkin 
equation ([23|) . as expected. Although not a martingale, properly modified, 
it converges to the density matrix predicted by Von Neumann measurement 
theory. 

The article is organized as follow : In section [2] we define the repeated 
QND measurement process we study. In section[3]we establish the link with a 
classical random process in which the pointer state distribution is repeatedly 
updated through Bayes' law. In section [J] we prove the convergence of the 
pointer state distribution under some assumptions and we determine the 
convergence rate in general cases. In section [5] we extend these results to the 
degenerate case. Finally in section [6] we study the time continuous scaling 
limit of our model. Some technical details appearing along the article are 
postponed to appendices. 

2 QND measurements as stochastic processes 

The aim of this section is to describe the relation between repeated non- 
demolition measurements, positive operator valued measurements (POVM's) 
and classical stochastic processes. 

2.1 Repeated indirect quantum measurements 

Let us consider a quantum system with initial density matrix p. Repeated 
non-demolition measurements aim at getting indirectly information on the 
system (without demolishing it as a projective quantum measurement a la 
Von Neumann or a direct connection to a macroscopic apparatus might do). 

To gain information, we let the system interact with another quantum 
system called the probe, and then perform a Von Neumann measurement on 
the probe. Assume the probe is initially in the pure state |^')(^|. Let U 
be the unitary operator, acting on the tensor product Hilbert space ^sys. 
"Wprobe, coding for the interaction between the system and the probe. After 
interaction, the system and the probe are entangled. Their joint state is 



A perfect non-degenerate projective measurement is then performed on the 
probe. That is, one is measuring an observable with a non-degenerate spec- 
trum i S X. Let {\i)} be the corresponding eigenbasis of "Hprobc- If the 
output of the probe measurement is i, the system state is projected into 
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because the probe and the system have been entangled. This projection 
occurs with probabihty 

tt{i) := Tr[{i\U\^)p{-^\U^i)]. 

We do not have to worry about cases where 7r(z) = 0, because these cases, 
almost surely, never happen. 

The process of "interaction plus probe measurement" is an example of a 
positive operator valued measurement (POVM). Let us define operators Mj, 
acting on the system Hilbert space, by 

Mi := {i\U\^). 

They satisfy M^Mi = Igys. as a consequence of the unitarit}^ of U . After 
measurement with output i, the system density matrix p'{i) can be written 
as 

p'{i) = ^MipM}, 

with 7r(i) = Tr[Mj pMj]. This characterizes a POVM. 

Let us now assume that we repeat the process of "interaction plus probe 
measurement" ad libitum. As we shall see below, even for purely practical 
reasons, it is useful to keep the freedom of changing some or all features of 
the process. For instance, the experimenter might tune (or let fluctuate ran- 
domly, or tune but leaving a certain amount of randomness or ...) the initial 
state 1^*) of the probe at each step. Or he/she might tune (or let fluctuate 
randomly, or ...) the interaction operator U at each step, for instance by 
playing on the time lapse that the probe spends close enough to the sys- 
tem to interact significantly with it. He/She might even tunc (or ...) the 
type of probe (in particular the dimension of its Hilbert space) at each step. 
Finally, he/she might tune (or ...) the non-degenerate probe measurement 
(equivalently the ^probc basis made of its eigenvectors 

We let l^'n), Uji, In denote the initial state, interaction operator and set 
of possible outcomes of the n*'' step. Setting po := Pi '■= p' Stud so on, we 
get a random recursion equation, namely that for i Ein- 

Pn = -^M,("W_iM,(")^ (1) 

with probability 7rn(i) = Tr[Mf"Vn-iMf"^^] where M^"^ = (i|;7„|^'„) (note 
that the meaning of the expectation itself, i.e. the Hilbert space with respect 
to which it is taken, might depend on n, however we may arrange to choose 

the X^'s so that i determines {i\ completely). 

It is worth noticing that, under such an evolution, a pure state remains 
a pure state, that is: ](/)„) = M-'-"'''|(/)„_i)/||M-'-"'|(/!).„_i)|| with probability 

n—i)- This case is included in that of density matrices. 
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Let us now specialize this scheme in such a way that it preserves a pre- 
ferred basis of the system Hilbert space. That is, we assume there exists a 
fixed basis {\a)} of T^sys. such that ah interactions can be decomposed as 



where the f7„(a)'s are unitary operators on "Hprobe^- The states |a) are cahed 
pointer states. The density matrices |Q)(a| with |q) a pointer state are fixed 
points of the recursion relation ([T]). 

The pointer states have to be eigenstates of the system Hamiltonian 
Hs for these indirect measurements to be quantum non-demolition (QND) 
measurements, since there is a waiting time between two successive indirect 
measurements during which the quantum system evolves freely. That is, 
Hs = -E'«|cK)(a| where is the energy of the pointer state \a) for the 
free system. 

After each indirect measurement one gains information on the system 
state. Repeating the process (infinitely) many times amounts (as we shall 
explain) to perform a measurement of a system observable whose eigenstates 
are the pointer states. This observable commutes with the system free evo- 
lution. A system in one of the pointer states remains unchanged by the 
successive indirect measurements. 

It has been shown in ref.|T] that a system subject to repeated QND 
measurements as described above converges toward one of the pointer states. 
This convergence was only proved in the case where the probes, interactions 
and observables on the probes are all the same. A non-degeneracy hypothesis 
was also used. One of the present article aim is to generalize the convergence 
statements without those assumptions. 

A word on terminology: we are going to name partial measurement one it- 
eration of "interaction plus probe measurement" and complete measurement 
an infinite sequence of successive partial measurements. 

2.2 A toy model 

We shall illustrate this framework and our results with a simple toy model 
inspired by experiments done on quantum electrodynamics in cavities [2]. 
The present work is actually inspired by these experiments. There, the 
system is a monochromatic photon field and the probes are modeled by two 
level systems. The observable we aim at measuring is the photon number. 
This is a non-demolishing measurement. 

The system-probe interaction is well described by the unitary operator 



where p is the photon number operator, e the energy of a photon and At the 
interaction duration. 




(2) 



a 



U = exp[ -i(eAtp 0h + -rp(S> 0-3) ] 
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This interaction amounts to the rotation of the two level system effective 
spin half if the cavity happened to be in a photon number operator eigenstate. 
The probes are assumed to be initially in a state 1^') = e"*'^'^3|+i) where 
l+i) is an eigenvector of ai corresponding to the eigenvalue +1. The probe 
observables, which are measured after the interaction between the system 
and the probe has taken place, are Oq' = e~*^ °"''(T2e*^ Their eigenvectors 
are \±'g) = e-^'^'^3|±2). 

The resulting POVM operators for the process of "interaction plus probe 
measurement" only depend on the difference between the two angles 9 and 
e'. They are M^"^' with 

_g-ieAt p ^_|_^ I gi(6»' -0)0-3 - if p(g)o-3 I _|_^ ^ 

Using the identity e"*^'^^ = cos(0)l2 — i sm{9)a^ and cr3|+i) = | — 1), one gets 
M^-^' = -L e-i(^Atp±f ) [ _ ^/ ^ !^p) ± sin(0 -9' + Ip)] 

One may verify that M^^^'Uf^"^' + M^"^''^M^"^' = D. Remark that if 
G N is a fixed photon number state, then |(p|M^|p)| is identical to 
\{p + 4k\M^\p + 4A;)|, with fc G IKI. This property leads to degeneracies in 
iterative QND measurement methods. Two states \p) and \p + 4/c) can not 
be distinguished by this method. These degeneracies are discussed in section 

EE 



3 Measurement apparatus and Bayes' law 

The aim of this section is to reformulate (part of) of the iterative QND mea- 
surement method in terms of classical probability theory. We are interested 
in the eigenstate probability distribution qn{oi), with 

qn{a) := {a\pn\a), 

and its evolution during the iterative procedure. At each step, the system 
density matrix is updated via the relation ([TJ, and as a consequence of the 
factorization property ([2]), this distribution is updated through the random 
recursion relation: 



qn{a) = qn-i{a)- 



|MW(i|a) 



E;3 9n-l(/3)|MW(i|/3)|2 

with probability X]/? 9n-i(/5)|M(")(i|/3)|2. We have defined M(")(i|a) 

{i\Unia)\-i'n). 
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Since |M(")(i| q;)P is the probability to get a probe measurement result 
i conditioned on the system state being \a), we introduce a (hopefully) sug- 
gestive notation 

Pn{i\a) := |M(")(i|a)|2 

We have Z]jPn(^|«) = 1 since X;iM(")(z|a)tM(")(i|a) = n. The recursion 
relation on the distribution reads 



with probability 7rn{i) = Qn-i{P)Pn{i\P)- This update rule corresponds 
to Bayes' law. The study of the eigenstate distribution convergence is thus 
a question of classical probability theory. 

Let us now put on stage the classical probability theory framework we 
shall be using. 

We imagine building a measurement apparatus which performs a se- 
quence of partial measurements. As we have stressed above, we allow for a 
protocol where the characteristics of partial measurements may vary at each 
step. However, we shall assume that these characteristics arc chosen within 
some finite set O called the set of measurement methods. In a quantum set- 
ting, one measurement method o is a triplet (probe state 1^*), interaction U, 
probe eigen-basis of direct measurement). Each measurement method 
o E O defines a set, called the spectrum of o and denoted by spec(o), of pos- 
sible outcomes. For each o we have a family of probability measures p°{-\a) 
on spec(o) indexed by a G 5, where S is the index set of pointer states. 

As time goes by, the experimenter records the sequence 01,11,02,12, ■ " 
where oi is the first measurement method, ii the outcome of the first partial 
measurement, performed using method oi, and so on. So it is natural to 
take as the space of events the space Q of infinite sequences (oi, ii, 02, ^2, • ' ' ) 
where each o„ belongs to O and each z„ to spec(o„). To be even more formal, 
set E := Uoeojo} x spec(o), so that E is the set of couples {o,i) with o G O 
and i £ spec(o). Then := E^ . 

For a finite sequence {oi,h,02, ■ ■ ■ ,in,On+i) e E"- xO, Bo^,i^,02,- ,in,Ou+i 
is defined as the subset of made of all those w's whose first 2n + 1 com- 
ponents arc oi,ii,02,--- ,i„,o„+i. Wc define Boj^^i^,02,- ,o„,i„ analogously. 
We let J>i be the cr-algcbra generated by all the -Boi,ti,02.- - .in-o„+i • Note 
that J^o is the a-algebra generated by all the Bg^, i.e. J-q codes for the 
first measurement method choice. For convenience we define T-i = {0,^2}. 
Then F := (J-"-!, J-q, J-i, • • • ) is an increasing sequence of cr-algcbras. We 
take to be the smallest cr- algebra on 17 containing all the mak- 
ing (yi,J-,f) a filtered measurable space. We could define another filtra- 
tion by taking 7"^ to be the cr-algebra generated by all the -Boi,ji,02,- -,on,i„- 
While this may seem superficially a more natural choice of filtration, we 
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shall see below that F is slightly more convenient. There is a natural col- 
lection of measurable functions on {il.,J^), namely the projections : for 
w = (oi,zi,02, • • • ,On,in,---) wc Set On(w) = o„, I„(a;) = in- These can 
be used to define counting functions that play an important role in the fol- 
lowing. We set e„(o,i) := lo„=oJ„=i and Nn{o,i) := T.i<m<n^rnio,i) (with 
the usual empty sum convention N(){o,i) := 0). 

The first task is to put a probability measure on Q. The next one will be 
to define a sequence of random variables on O solving the recursion relation 

If the measurement methods are given, the distributions of partial mea- 
surements are described by the p°{-\a). So what remains to be discussed is 
how the measurement methods are chosen, and we put the condition that this 
does not involve precognition. We suppose that a collection of non-negative 
numbers = 1, do(oi) (for all oi G O), • • • , dn{oi,ii, 02, • • • ,i„,o„,+i) (for 
all oi € O, ii G spec(oi), • ■ • , On+i G O) is given in such a way that 

c?„(oi,ii,02, • • • ,i„,o„+i) = d„_i(oi,ii,02, • • • ,Zn_l,0„). (4) 

o„+ieO 

It is not difficult to produce such a collection. For instance if Co(-) is 
a probability measure on O and for each n > 1 c„(-|o, i) is a probability 
measure on O indexed by (o, i) € E then 

Cl (02 |oi , il) ' ' ' Cn {On+l\On, '•n J 

does the job. We call this special choice the Markovian feedback protocol. 
Some special cases are of interest. If we assume that for n > 1 c„(-|o, i) = 
Cn(-) does not depend on o,i, we arrive at something we call the random 
protocol. On the other hand, if 60 is an element of O and 5„, n > 1 a family 
of maps from E to O, then taking Co(-) := S.^bo and, for n > 1, c„(-|o, i) : = 
^■,bn{o,i) arrive at the description of an experimenter deciding of the next 
measurement method by taking into account the previous measure outcome. 

Given a collection of such non negative numbers (iy^(o]^, * * * ; ^n? ^n+i)? 
and using the Kolmogorov extension theorem, it is easy to see that there is 
a unique probability measure Pq, on {^,J-) such that 

Pa(^oi,ji,02,-,j„,o„+i) = • ■ ■P°"{in\a)dn{0i,ii,02, ■ ■ ■ ,i„,0„+i). 

Indeed, we see that the mandatory consistency condition is fulfilled : if the 
left-hand side is summed over o„+i (using (jj])) and then over i^i the forniulci 
for PaiBoi,ii,02,-,in-i,on) is recovered, which is needed since _Boi,ii,02,- ,j„-i,o„ 
is the disjoint union of the -Boi, n, 02, ■■■ ,«n, 0,1+1 '^^^^ possible On+i and i„. 
The normalization condition Pa{^) = 1 and the positivity condition are 
obvious. 
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Note that in general, conditional on the sequence of measurement meth- 
ods oi,--- ,o„, one has Pq,(?i,--- ,in\oi,--- ,o„) / p°i(ii|a) • • •p°"(i„|a). 
This is due to the feedback. For the cases when the d„'s do not depend on 
the outcomes, in particular for the independent random protocol, equality is 
recovered. 

We define also, 

P=Y,Qo{a)Pa. (5) 

We use Eq and E to denote expectations with respect to Pq and P respec- 
tively. 

A simple computation shows that, for each a, the conditional probability 

p (n - r, \n T n T^- (^n{0l,h,02,- ■ ■ Jn,On+l) 

^aV^n+l — On+l\Ui, li, • • • , In 



dn-l(Oi, /i, O2, • • • ,/„_!, On) 

whenever the denominator in nonzero. The same formula holds for P. The 
right-hand side is simply c„,(o„+i|On, I„) for the Markovian feedback pro- 
tocol. So indeed, the functions do, di, • • • embody the probabilistic de- 
scription of the choice of measurement methods. 

These definitions may seem arbitrary at that point, but now we can make 
contact with the initial problem. Define a sequence of random variables 
Qn{o<) by the initial condition Q{){a) = qo{a) and the recursion relation 

^ f ^ Qn-i{a)p^"{In\a) 
Qn{a) - 



To show that the recursion relation Q is verified, we need to show that the 
transition probabilities are correct. A simple way to do that is to solve this 
random recursion relation. For uj S -Boi,ji,02,--- ,jn,o„+i one checks that 

^ ( ^ _ 9o(«)p°' {h\a) ■ ■■p°"{in\a) 



whenever P(i3oi,ji,02,--- ,«n,o„+i) 7^ 0. Note that this condition ensures that 
the denominator of Qn{ct,uj) is nonzero. We observe that whenever defined. 
Quia, id) > and Yl^Qnia,uj) = 1. If PiBo^^i^,02,- ,i„,o„+i) = the value 
of Qniot,uj) is mostly immaterial from a probabilistic viewpoint, because in 
any case the full sequence Qn{ct, is well-defined on a set of P-measure 1 
(note that the collection of all Bo^^ii^o2,--- ,i„,o„+i is countable, so the collection 
of those with P-measure is countable as well, or empty). 

Since there is no dependence on o„+i on the right-hand side we observe 
that for a; G 0, and conditional on J-n-i, 

Qn{a) - 
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with probability 



So the recursion relation ([3]) is recovered with the identifications g„(a) — )• 
Qn(cK) and — ^ To summarize, we have accomplished our goal : find 
a probability space on which ([3]) has a solution, which we have even written 
explicitly. 

In the following sections we study the convergence of Qn{-) and its de- 
pendence with respect to the initial pointer state distribution. On our way, 
we shall understand the probabilistic meaning of the recursion relation ([3]). 



4 Convergence 

In [1] the convergence of (?n(-) has been shown under the hypothesis that 
only one partial measurement method o is used. The properties of the 
limit were elucidated under the further assumption that for every couple 
of pointer states (a, /3) there exists at least one partial measurement result 
i such that p°{i\a) ^ p°(i|/3). This last assumption can be understood as a 
non-degeneracy hypothesis because two different pointer states a, /3 do not 
induce identical partial measurement results distribution p°{i\a). Our aim 
is to generalize the convergence of Qn{-) while weakening the hypotheses 
made in pj . We discuss the convergence when different partial measurement 
methods are used. We focus on the influence of this extension on the rate 
of convergence. The degenerate case will be studied in section 15.11 In the 
case of one measurement method, a convergence result similar to that of 
[T] has been obtained by H. Amini, P. Rouchon and M. Mirrahimi through 
sub-martingale convergence in 



4.1 Convergence with different partial measurement meth- 
ods 

The extension of the convergence result of [T] to cases with different mea- 
surement methods is straightforward. From the fact that, conditional on 

•f' n—li 

f) ( )^ Qn-l{a)p"" (in|«) 

with probability '^i^^g Qn-iW)p°"i'''n\/3), the average of Qn(o)) again con- 
ditioned on J-fi^i, is 

E Qn-l(/3)p°"(in|/3) = Qn-l(a). 
i„especo„ 13 
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So Qn-i{a) is conserved in average. Though the computation involved to 
prove it is essentiaUy the same, a mathematically cleaner statement is that 

HQn{a)\J^n~l) = Qn-l(a), 

i.e. each Qn{oi) is an F- martingale. 

In fact, Qn{ct) has a deeper probabilistic meaning, which makes the mar- 
tingale property obvious. 

For a while, forget the previous definition of Qn{ct). Observe that, under 
the assumption that qo{oi) > for every a £ S, any set of P-measure has 
also Pc-measure 0. Then the Radon-Nikodym theorem states that for each 
a £ S, there is a P-integrable non-negative random variable Q{a) on (0, J^) 
such that 

qo{a)Ea{X) = HQ{a)X) 

for every Pa-integrable random variable X on {0,,J^). The random variable 
Q{a) is a Radon-Nikodym derivative of qQ{a)Pa with respect to P. It is 
obvious that two Radon-Nikodym derivatives can differ only on a set of P- 
measure 0: in that sense the Radon-Nikodym derivative is unique if it exists. 
We have also that, P-almost surely, Q(a) = 1, so that, P-almost surely, 
each Q{a) < 1. This existence theorem is a bit abstract but if one replaces 
J- by J-'n one can get a concrete formula. The same argument ensures the 
existence of a P-integrable non-negative random variable Qn{c() on (Q,J^n) 
such that qo{a)Ea(X) = E{Qn{a)X) for every Po-integrable random variable 
X on (r2,J-'„). As J-'n is finite, it suffices to let X run over the indicator 
functions for the -Boi,n,02,--- ,in,o„+i • This implies that 

Qn{a,uj) = qo{a)- 



for every uj E Bo^^ii^o2,--- ,i„,o„+i (such that the denominator is nonzero, else 
the value of Qn(cy,U!) is immaterial). Explicitly one finds 

qo{a)p°^ {ii\a) ■ ■ ■ p°'' {in\a) 



E^esQoil3)p°^iii\f3)---p°"iin\/3) 

on -Boi,«i,02,--- ,i„,o„+i • 

This is exactly our previous definition of Qnict), which is probabilistically 
a Radon-Nikodym derivative. This makes the martingale property obvious 
without any computation, just by the definition and general properties of 
conditional expectations. In fact, Qn{c() is a closed martingale: 

Qn{a) = E{Q{a)\Tn). 

As Qn{ct) is also bounded, the martingale convergence theorem [5] ensures 
that Qnict) — )• Q{a) P-almost surely and in L"'^(17, T, P). Since by assumption 
5 is a finite set, all Qn{ctys converge simultaneously P-almost surely. 
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Let us now observe that the fohowing two statements are equivalent: 

— the measures Pq are ah mutually singular, 

— there is a collection {^}a)a£S of disjoint measurable subsets of Q such 
that, P-almost surely 1^^ = Q{a). 

The proof is simple. The statement that the measures are all mu- 
tually singular is equivalent to the existence of a collection {fla)aes of dis- 
joint measurable subsets of il. such that P^(r2Q) = (5q,^^ for all a,/3 S S. 
From the defining property of Radon-Nikodym derivatives, lfi^Q{a) is also a 
Radon-Nikodym derivative of go(Q^)Pa with respect to P, and E{Q{a)ln^) = 
qQ{ct)Sa,i3 which by positivity implies that, for a ^ f3, (2(a)lr2^ = except 
maybe on a set of P-measure 0. Hence P-almost surely l^^Q{f3) = Q{(i)Sia^p. 
Summing over (3 gives In^ = P-almost surely. The converse is also true: 

if there is a collection {^a)a&s of disjoint measurable subsets of such that 
P-almost surely 1^^ = Q{a), then the measures P^ are all mutually singular 
and concentrated on the Q^s. 

A striking consequence is that, if the measures P^ are all mutually sin- 
gular, for each u; in a set of P-measure 1, Qn(cK) converges to either or 1, 
and it converges to 1 with probability P(Oa) = qQ{a). 

Hence when the measures Pq, are all mutually singular there is a full ex- 
perimental equivalence between an infinite sequence of partial measurements 
and a direct projective measurement on the system. We further study this 
equivalence in section W% 

We shall now give a criterion, that the experimenter may enforce on the 
protocol, ensuring that the measures P^ are all mutually singular. This 
involves a non-degeneracy hypothesis, similar but weaker than that made in 

m- 

We say that o G O is recurrent in a; G if 0„(a;) = o for infinitely many 
n's. Our (sufficient) criterion for all the measures P„ to be mutually singular 
is: 

There is a subset Os of O such that 

- Each o G Os is recurrent with probability 1 under each Pq, 

- For every a,/3 G 5, a 7^ /3 there is some o G and i G spec(o) such 
that p°{i\a) /p°(i|/3). 

This condition says that with probability one, infinitely many partial 
measurements that distinguish between any two states of the system will 
occur. 

Consider the event Aq made of the w's such that Qn{oi) converges for 
each a and o is recurrent. Note that by our assumptions P(y4o) = 1- We 
show that for any i G spec(o) 

Q{a)Y,Q{l)P°m = QW{i\a) 
7 
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on Ao- There are two cases to consider. Either {On{^^), In{^)) = {o,i) for 
infinitely many n's: then the announced relation follows by taking the limit of 
the basic recursion relation along a subsequence. Or (0„(cl!), /^(a;)) = (o, i) 
for only finitely many n's: then, as shown in Appendix lAl Q{l)p°{Al) — 
so that in particular Q[a)p°{i\a) = and the announced relation still 
holds. 

This implies that 

Va,/3 G 5, Q{a)Q{l3){p° {i\a) -p°{i\l3)) = on A„ for every i e spec(o). 

(7) 

Then, by ([7]), Va, /3 G 5, a / /? one has Q{a)Q{f3) = on Aq^ := Cio^Os^o, 
which has P-measure 1. As the sum of the (5(a)'s is 1, this means that on 
Aq^ one has Q{a) = 5a,r where T(uj) is some w-dependent element of 5. So 
there is a family of disjoint subsets f^Q, of Ao^ such that Da^a has P-measure 
1, and Q{a) = except maybe on a set of P-measure 0. 

We shall give two examples. 

For the first one, the task to ensure that enough measurement methods 
o are recurrent is left to the experimenter. 

The second example is the Markovian feedback protocol. For this pro- 
tocol, under Pq,, the process (o^^, i^^) is Markovian with transition kernel 
Ka,n{o,i',o' = p° {i'\a)cn{o'\o^i), with initial distribution p°{i\a)cQ{o). 
Recurrence questions are well under control at least when the kernels do 
not depend on n. So we assume that c„(o'|o, i) = c{o'\o,i) is time in- 
dependeniH and set Ka{o,i;o' = p°' {i'\a)c{o'\o,i). The product struc- 
ture of Ka leads to introduce the reduced transition kernel K^'^{o]o') : = 
X]iespec(o) 0- '^^^ ^^^y classical Markov chain computa- 

tions to make sure that the measures Pq, are all mutually singular. Assuming 
that the reduced Markov chain K^'^'^ is irreducible and aperiodic, it admits 
a unique invariant probability fi^'^ on O, which is strictly positive. Then 
all partial measurement methods with ^fa'^{o) > will be recurrent on a 
set of P^-measure 1. Moreover the full Markov chain has a unique invariant 
probability ^a{o,i) = p° {i\a) ^^^'^ (o) . Then the strong law of large numbers 
for Markov chains states that Nn{o,i), the number of occurrence of {o,i) up 
to the n*^ experiment, satisfies 

Nn{o,i) ~ n^a{o,i) for large n 

on a set of Pa-measure 1. This ergodicity result will be put in use in section 

To summarize, we have proved that if there are enough P^ recurrent 
partial measurement methods then the measures P^ are all mutually singular, 

^Were c„ periodic in n we could look at a Markov chain with a larger state space to 
reduce to that case. 
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so that there is a fuh experimental equivalence between an infinite sequence 
of partial measurements and a direct projective measurement on the system. 

4.2 Conditioning or projecting 

In the previous section we pointed out the connection between complete 
measurements and direct projective measurements. This holds whenever the 
measures Pq, are mutually singular. 

Under this hypothesis, we show in this section that we can solve the 
random recursion relation jS]) on a space where the final pointer state is de- 
termined before the measurement process starts: for the class of experiments 
we are dealing with, it is consistent to assess that the total measurement out- 
come can be decided in advance and by a classical probabilistic choice. 

We assume that the measures Pq, are mutually singular, and, as usual, 
that all qo{a) are > 0. To avoid clumsy statements, we remove from 17 the 
(P-negligible) set of events for which either (5n(a) is not defined for all n's, 
or Qnict) does not converge to or 1. So we assume that the sets Qq. form 
a partition of 0, and the random variable T defined by lim„_j.oo Qn(ct,a;) = 
'^ujeUa — ^a,T{uj) is defined everywhere on 0. We let A be the smallest a- 
algebra making any measurable. We claim that if X S L^{Q,T, P) is any 
integrable random variable, 

E{X\A) = Et(X). 

This can be rephrased as : conditioning P on the limit of Qn{ct) being 1 
leads to Pq,. This is essentially obvious from the Radon-Nikodym viewpoint 
and the fact that the measures Pq, are mutually singular. But a direct com- 
putation is easy. The fact that S is finite (countable would do the job as 
well) has two consequences. First, any ^-measurable random variable Y can 
be written as a linear combination Y = ^Qg^yalrJa- Second, to test that 
E{X\A) = Y it suffices to check that E{XlnJ = E(YlnJ for every a. Now, 
by definition, E{XlnJ = qo{a)EaiX), whereas E{YlnJ = qo{a)ya, so, if 
E{X\A) = Y,y^ = Eq(X), i.e. Y = Eae5 = Et(X). Hence 

E(X|^) = E-x{X) as announced. 

This proves the equivalence between projecting first on a given state a 
and conditioning on the limit state being a. 

4.3 Convergence rates and trial distribution independence 

Experimentally, the initial distribution go(") may not be known. One would 
then use the sequence of partial measurements to gain information and re- 
construct it from these measurements. This may be done using Bayes' law 
starting from a trial distribution go(") (supposed to be nowhere vanishing) 
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and recursively improving it using the relation 



P°"(in|a) 



Zl3qn-lil3)p°-iin\ 



if the outcome is i, which happens with probability g„_i(/3)p°" 
The difference with eq.Q is that the recursion involves ^n(') and not qni')- 
However, the probability is the one given by the (?„(•). If the initial trial 
distribution ^o(") coincides with the initial system distribution qo{-), then 
Qn{') = Qni') for all n. Both gn(') and qn{-) are realization dependent. We 
shall define the random process Qn{') as Qni') in dS]) but with a different 
initial distribution 

Qn{a) ■■= go(a); 



E/3^o(/3)p''Kn|/5)---p°"(in|/3) 

The probability law still depends on the true initial distribution. Notice that 
Qni') is not an F-martingale under this law, contrary to Qni')- As we shall 
show, they nevertheless have identical limit, that is: lim„_i.oo Qn(a) exists 
and is equal to 1^^ with P-probability 1. 

Moreover, if a time independent Markovian feedback protocol is used, 
the convergence of the state probability distribution is exponential. Its con- 
vergence rate is the mean relative entropy of the partial measurement result 
distribution conditioned on the system be in the state T with respect to the 
one conditioned on the system be in the state a. This means that for n large 
enough, 



with 



Qnia) ~ e-"^(^l"), for a / T (8) 

p°ii\(3 



Sif3\a)'.= Y,f^T\o) P"i'\f^)^^ 

oGO i6spec(o) 



p°ii\a) 



(9) 



Here, all the p°ii\a) are assumed to be strictly positive, thus any (o, i) with 
o £ Os is recurrent. In the case of time independent random protocols, the 
rate is the same with fi^'^i') replaced by c(-) the distribution of measurement 
methods. This coincides with the result of [Ij if O contains only one partial 
measurement method. 

The independence of the limiting distribution with respect to the initial 
trial distribution is obtained whenever one starts with a trial distribution 
such that qoia) > wherever qoia) > 0. This happens for example if we 
start with qoia) > for any a £ S. 

To see this, we analyse the behavior of Q„ under the probability measure 
P := X]ae5 9o(ck)P«) which can be seen as a trial probability measure on Q 
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as it corresponds to a system initialy in the trial state. Under P, is a 

martingale, so by the above arguments, it converges P almost surely to 
As by hypothesis qo{a) > whenever qo{a) > (which can be rephrased as: 
P is absolutely continuous with respect to P), a subset of $7 of P probability 
1 has also P probability 1. So lim„^oo Qn(o) = If^^ = hm„_j.oo Qn{ct) with 
P probability 1. 

What is less direct is the determination of the convergence rate. This 
requires controlling the behavior of the counting processes Nn{o,i). As re- 
called at the end of section l4.1l Nn{o, i)/n ^ fJ-a{o, i) on a sel|f|of Pa-measure 
1. We want to infer that P-almost surely, 

lim Nn{o,i)/n = fir{o,i)- 
To prove it, we set C{o,i) := {lim„_j.oo A'„(o, i)/n = /Ut(o, i)} and write 

a a 

where we used 



P(£(o,i)|J7„) = P( lim Nnio,i)/n = fiUo,i)\^a) 

n— >oo 

= Pa{ lim Nn{o, i)/n = ^^(o, i)) = 1. 

n— ^oo 

Observe now that Qn{o) can be expressed as a function of the counting 
processes: 



Qn{oi) 



E/3^o(/3)n(o,)ei^P°(^l/3)^"(°'^) 



Under the hypothesis that all p°(i|a)'s are >0, we have that Qn{o) > for 
every n, and the logarithm of the ratio between Qnict) and Qn{P) is well 
defined. Using the previous result, we have 



lim — In 

ra— >oo n 



QnW) 

Quia) 



p"{i\a) 



J2 MT(o,i)ln 

{o,i)&E 



p°(z|T) 



+ In 



P°(^|T) 
p"{i\a) 



^ g-n5(T|/3)gn5(T|a) 



Then, for a large enough n, 

Qn{f5) 

Qn{a) 

^For the record, the set {N„(o,i)/n — >■ fia{o,i)} is measurable, as it can be written 
n^eN* UnoSN n„>„o {'^ £ |A^n(o,i;a;)/7i~p„(o,i)| < 1/m}. 
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with S{T\a) the mean relative entropy, S{T\a) = X^o^o A*t'^(o)S'°(T|q;) 
where 

S%T\a)= Yl P°(^|T)(ln[p°(^|T)]-lnb°(^|a)]) 

iSspec(o) 

The relative entropy is always non negative, subsequently, the mean relative 
entropy is non negative too. Moreover, the mean relative entropy is null if 
and only if T = /3 (all relative entropies null). 

Using this property and QniP) = 1, we obtain for a 7^ T, 

Then, to leading exponential order 

Hence, for n large enough, we proved that 

^ r 1 _ if a = T 

"^"^"^ - 1 const.e-"^(^l°) else 

The limit distribution docs not depend on the trial initial distribution but 
only on the complete measurement realization. The probability to have 
Qoo{oi) = Sa,^ equals qoil)- With time independent Markovian feedback 
protocol, the convergence is exponential with a leading rate S{T\a). 



4.4 Convergence rate tuning 

Most of the time, when performing a measurement, one prefers it to take 
as little time as possible. The use of different partial measurement methods 
allows us to tune the convergence rate. Let us take an example. Suppose 
we want to discriminate between three possible pointer states of a system, 
and suppose that the partial measurement methods give only True/False 
as possible outputs. We denote T, F the partial measurement results and 
1,2,3 the pointer states. Each partial measurement method can be tuned 
to maximize, up to a measurement error e ^ 1, the probability of one of its 
outcome knowing the system is in one of the three pointer states. We shall 
show that this is not enough to maximize all convergence rates for arbitrary 
limit pointer state. It is the use of different measurement methods picked 
randomly that allows us to overcome this convergence rate problem. 

Let us consider for instance two measurement methods. The first one, 
denoted a, has conditioned probabilities 

p«(r|l) = £,p"(r|2) = g,p«(r|3) = l-e 
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with q = 0(1). The second one, denoted b, is obtained by switching the 
probabihty conditioned on 1 and 2, that is 



p\T\l) = q,p\T\2) = s,p\T\3) = l-s 

Let us now look at the convergence rate conditioned on the hmit pointer state 
to be 1. These are coded in the relative entropies. If only the measurement 
method a is used, one has: 



S"(l|2) =eln 
5"(1|3) =eln 



+ (1 - e) In 



1 



1 



+ (!-£) In 



l-e 



~ -ln[l -g] = 0(1) 
- ln[e] » 1 



The convergence of Qn(3) toward is quick but the one of (5n(2) is rather 
slow. If measurement method b is used the interesting relative entropies are 
now 



5''(1|2) =gln[^J +(l-g)ln 



1 



S\l\3) =qln 



l-£ 



+ In 



^ — ln[s] 



ln[s] 



All the convergences rates are then high if the limit pointer state is 1. But 
if the limit pointer state is not 1 but 2 then, using only the measurement 
method b, the relative entropy 5^(21 1) is 



5^(2|1) =£ln 



(1 - e) In 



l-e 



1-g 



Hq] = 0(1) 



and the convergence rate toward 2 is slow. 

Now, if at each time one of the two measurement methods is used with 
equal probability ^. The convergence rate for any i,j with i j is 

sm = ^is^m+s'm) ^-hs] » i 

As a consequence, the convergence rate is always high, whichever the limit 
pointer state is. 

In the toy model, if the first partial measurement method correspond to 
9-9' = ^, then 

5^(013) ~ 0.116 

This is the slowest of all convergence rates. If the partial measurement 
method with ^ — ^' = | is introduced and the partial measurement methods 
are chosen with equal probabilities each time, then 



5(0|3) ~ 1.18 
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and the slowest of all convergence rates is 

S{1\3) ~ 1.10 

If only the ^ measure is used and the limit pointer state is 0, then a the- 
oretical 99% confidence level is reached after about 50 measures. With the 
use of the two different partial measurement methods the same confidence 
level for the same limit state is reached in 5 measures. The same number of 
measurements is needed if the limit pointer state is 1, 2 or 3. 

5 Degeneracy and limit quantum state 

Often the quantity we measure is a property common to several pointer 
states. In the quantum case, this corresponds to a degenerate projective 
Von Neumann measurement. There, at least two different eigenstates share 
the same eigenvalue. For our measurement process, degeneracies happen 
when several distributions p°{-\-) are equal for different pointer states, so 
that some states cannot be distinguished. For example, in our toy model, 
whatever 6 — 9' is, we have p^~^ = p^~^ + with k an integer. 

The pointer state with p photons cannot be distinguished from the one with 
p + 4k photons. 

In this section we study the system state evolution when degenerate 
repeated partial measurements are performed. In a first part we show the 
convergence of the system pointer state distribution. In a second part we 
focus on the quantum case and the influence of phases introduced between 
degenerate states by the repeated partial measurement process. 

We shall partition the set of configurations into sectors. Let us define 
an equivalence relation among pointers by identifying two pointers whose 
partial measurement distributions are identical. That is: two pointers a and 
/3 are said to be equivalent (denoted a ~ ^) if, for any partial measurement 
method o and result i, 

p%i\a)=p%i\/3). 

By definition the sector a is the equivalence class of a. In the toy model the 
sectors are the sets p = {p + 4A;, A; G N} with p = 0, 1, 2, 3. 

5.1 State distribution convergence 

We first look at the convergence of the pointer state distribution Qni') in 
case of degeneracy. The system distributions Qn(') induce probability dis- 
tributions Qni') on sectors by 

Quia) := J2 Qnia'). 
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Since sectors form a partition of the set of pointer states, we have Qn(a) = 
1. The initial probabihty of a sector is qoia) = ^^a'&a^oi^')- '^^^ recur- 
sion relation ([3]) can obviously be lifted to a recursion relation for the sector 
distributions, 

It is identical in structure to eq.(l3]) but with the bonus that it now is non 
degenerate. Lets assume that for two different sectors, it exists at least 
one P-recurrent partial measurement method distinguishing between the two 
sectors: if a 9^ /3, it exists o £ Og and i G spec(o) such that p°{i\a) 7^ p°{i\l3). 
Thus we can use the non-degenerate case results but applied to the sector 
distribution. Hence, Qn{') almost surely converge and 

Q{a) = Sa,r 

with T the realization dependent limit sector. The probability that the limit 
sector be 7 is equal to 90(7)- 

From the martingale property, the state distribution converge (not only 
the sector distribution), and the point which remains to be discussed is what 
is this limit. Thanks to the relation qo{a') Qn{ct) = <?o(q^) Qn(ck') valid for 
any n if a ~ a', we shall show that this limit is 

= 1 .o(«)/.-o(T) ifaGT (10) 

with T the limit sector. 

Indeed, the state distribution satisfies the recursion relation ([3]). Thus, 
if (^fc)fc=o,--- ,n-i cire the n first partial measurements results, one has 

Qn(a') _ m=iP'''(.ik\a) 



for any a' in the sector a. The right hand side only depends on the sector 

Qn(a) ^ Qn{a') 
qo{a) qoia') 



a, and thus '^"/".^ = if a' G a. From this equality, it follows that 



Q{a') = Q{a)- 



qo{a) 

Since, Q{a) = 6a,r, we have Q{a) = if q ^ T and 1 = Yl^'eT Q('^') 
90(0) 



■^feyQ'o(X), for any a G T. Hence, 



go(«) 

Q[aj = —^, foraGT 

The probability of convergence to a sector as well as the limit state 
distribution (|10p coincide, in quantum mechanics, with what would have been 
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predicted by Von Neumann rules for degenerate projective measurements. 
The approach we have been following so far, based on tools from classical 
probability theory, gives no information on the convergence of the density 
matrix off-diagonal elements. It is the next section's purpose to discuss 
the evolution of the density matrix p„ and not only the evolution of the 
probabilities (5n(a) = (a|Pn|«)- 

5.2 Density matrix convergence 

We are now interested in the convergence of the system density matrix. In 
the basis of pointer states we may write: 



with An{a,a) = Qn{oi). It evolves according to the recursion relation ([T]). 
The processes An{a^l3)^ are not martingales, and their convergence can not 
be obtain through the martingale convergence theorem. Actually, they do 
not always converge. To obtain convergence a unitary evolution process has 
to be subtracted. 

For each POVM, a phase between the pointer states is introduced by 
the operators M^°\ Even inside a sector this phase can be nonzero. This 
possibility comes from the degeneracy criteria we unraveled previously. Two 
pointer states, a, f3 can have a nonzero limit probability if they are in the 
same sector : a ~ /3. This criterion implies a norm equality |M(°)(i|a)| = 
|M(°)(i|/3)| for any i in the spectrum of any partial measurement, but not a 
fuh equality. So M(°\i\a) and M^°\i\j3) can differ by a phase. The density 
matrix converges either if this phase can be set to zero or if we absorb it 
through a transformation of the evolution. 

Let us write the operators M^°\i\a) in a phase times norm form 



Hp + Ha )\ct){a\. This is the most general Hamiltonian if one want U to 
fulfill the non demolition condition ([2]). In the above formula. At is the 
interaction time between the probe and the system. 

Let us define a unitary operator process, diagonal in the pointer state 
basis. 



Pn = Y,Ma,f3)\a){l3\ 



a,l3 





0(°)(i|a)Af„(o,i)) 



a) {a 



a 



and the unitary equivalent conjugate density matrix process 



Pn = 



UlpnU, 



n 
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The diagonal elements of p„, in the basis are not affected by this 

transformation. Their limits stay the same. Thus if a. or /3 are not in the 
limit sector T, according to the Cauchy- Schwartz theorem, Aoo(a,/3) = 0. 
We are then interested in the limit of the elements j4„(a, /3) with a,/3 G T. 

If /? G a, then go(o)^ra(tt) /3) = ao(cK) /?)Qn(o)- Repeating the discussion 
made in the section 15. 1| we get 



^oo(a,/3) 



ao(a,/3)/go(X) ifa,^eT 
else 



Hence, 'pn has an almost sure limit which coincides with the result of 
a Von Neumann measurement: poo is equal to pQ projected on the system 
subspace corresponding to the sector T. 

1 

lim pn = _ .^x Pt/OqPt (12) 

n-5-oo ~ ~ 

where Px := S^gx I7)('TI projector on the subspace corresponding to 

the sector T. 

In some cases the unitary operator process can be reduced to a deter- 
ministic one. For example if d^°\i\a) = e'-°\i\(3) for every (o, i) £ E and 
every a ~ /3, then we can chose C/„ = ^-i-^^Hs _ xhis is the case in the toy 
model, from the periodicity of trigonometric functions, in a given sector all 
phases introduced by partial measurements are equal. If the limit sector is 
p = {p + 4:k, k G IKI}, then the photon field limit state will be 

lim e-^*P/,„ e— = V ^o{p + P + Ak') ^ ^ 

One other example corresponds to the case where the phases 6^''\i\a) do 
not depend on (o, i) but depend on a. Then one can define -ffe//. = Hg + 
0{a)\a){a\ and Un = e-'^'^^^-ff- . 

In most of the cases the unitary evolution C/„ is a stochastic process and 
then in the limit n — t- oo, it remains a stochastic rotation inside the limit 
sector. When Un is deterministic, the remaining rotation is deterministic 
too. 



6 Continuous diffusive limit 

We shall now prove the convergence of the discrete processes we consider 
toward processes driven by time continuous Belavkin diffusive equations. 
Our proof, different from that used in |21) , allows us to derive the continuous 
measurement diffusive equation not only for the quantum repeated indirect 
measurement process but also for the macroscopic Bayesian apparatus we 
defined. The quantum case is a peculiar realization of it. 
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The time continuous equation is found as a scaling limit of the discrete 
evolution when n goes to infinity with t = n5 fixed [5 = Ai). We first study 
the pointer state distribution scaling limit, Qtict) := lim5_>o Q[t/(5] (ct)- The 
evolution equation for Qt{a) is given in eq. (fT5|) below. We then look at the 
time continuous limit of the density matrix evolution and get the Belavkin 
diffusive equation eq. ()23p . In the quantum case, the continuous limit requires 
rescaling appropriately the system-probe interaction Hamiltonian as 

H = Hs'S!)ip + h(^Hp + -^Hi (13) 

We present in some details the case with a unique partial measurement 
method. The results are then easily extended to cases with different mea- 
surement methods. 



6.1 Continuous time limit of the pointer state distribution 

We are here interested in the state distribution continuous time limit. The 
results presented in this section apply to the general Bayesian recursion rela- 
tion ([3]) — which in particular includes the case of repeated QND measure- 
ments. To begin with, we assume that there is only one partial measurement 
methocfl. Henceforth we suppress o from all the notations and let T stand for 
the index set of outcomes. Note that the two filtrations J-n and T'^ coincide 
and carry the information on the first n partial measurements. 

We assume that the conditional probabilities p{i\a) depend on a further 
small parameter 6, and are of the form 

p{i\a)=pa{i){l + ^5Ts{i\a)) (14) 

with pQ{i) > for all i's and that r(z|a) := lim^^o+ r5(z|a) exists. Then 
^^Po(^) = 1, so that the po(^)'s specify a probability measure, and for every 
(5, X^jPo (0^5(^1 a) = 0. The important point is that po{i) is independent of 
a. 

These hypothesis are of course satisfied in the quantum case with QND 
interaction Hamiltonian Hi = |«)(«| <S5 and rescaling Hj — )• -^Hj. 
Then po{i) = is the probability measure in absence of interaction, 

and 

We assume that ^ for all i. 

We first need to make precise the sense in which a limit on Qn{c() is to 
be taken. 

^This condition will be relaxed in section [6.1.41 
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For a fixed 00 ^ Q the limit lim5_>>o Q[t/S] (o^) is not expected to exisiH. But 
there is some hope that, properly defined, a limit for the law of the process 
Q[t/s]{c()j t G R"*" exists. We shall prove in Appendix iBl that this is the case, 
specifying a bit the kind of convergence that is involved. Under the limiting 
law, the process Qt{-) satisfies the stochastic equation 

dQtia) = Qt{a) ^ {r{i\a) - {ri)t) dXt{i) (15) 

i 

where 

Here Xt{i), with X^j^f(i) = 0, are continuous martingales with quadratic 
covariation 

dXt{^dXt{j) = dt {SijPo{i}-po{i)po{j)). (16) 

We shall show that a vector solving this equation is a bounded martingale, 
to which the martingale convergence theorem can be applied with results 
similar to those in the discrete case: 



Qoo(a) 



if a T 

qo{a)/qo{a) if a G T 



with T the limit sector. However, the sector definition is not the same as 
in the discrete case. In time continuous, a and /3 are in the same sector 
if and only if r(i|a) = T(i\/3) for all partial measurement result i. The 
probability for the system to be in the sector a in the limit t goes to infinity 

is qo{a) = Ea'ea9o(a')- 

The convergence is still exponential 

(5t(a) = exp (-t/rxa), if a T 
with characteristic convergence time t^jq^j 

= ^Po(i)(r(i|a) -r(z|7))' (17) 

i 

This coincides with the convergence rate we would have found by taking 
the relative entropy S{'y\a) scaling limit. However, it is somewhat difficult 
to decipher that it originates from a relative entropy by only knowing its 
expression in the continuous-time limit. 



Think of the simple random wallt: the convergence to Brownian motion is not sample 
by sample because has no reason to be close to 
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6.1.1 Preparation 

We work with the model {^,J-, P). 

Our derivation is based on the use of the counting processes Nn{i). Recall 
that No{i) = and that iV„(i) := Z^i<m<n ^m(^) n > 1, where : = 

is 1 if the n*^ partial measurement outcome is i and otherwise. 

We start by listing some properties of these counting processes and their 
relationship to the solution of ([3|). Then we shall formulate and prove the 
analogous statements for the continuous time limit. 

It is obvious that the filtration J^q,J^i,--- is the natural filtration of the 
vector counting processes A^^^. 

Also recall that the random recursion relation ([3]) can be solved in terms 
of the counting processes as 




A trivial but crucial observation is that under each P^, A^^^ is the sum of 
independent identically distributed (i.i.d) random vectors. 
As a first consequence, a simple computation leads to 

E j^gEtiE,A,(i)(JV„,(i)-7v„,_,(i))^ ^ 

k / 
a 1=1 \ i 

for A; > 1, arbitrary non-decreasing sequences of integers = no < ni < • • • < 
Uk of length k, and arbitrary (complex) A;(z)'s. A second consequence is that 
under Pq, each Nn{i) is a sub-martingale and Nn{i) = {Nn{i) — np{i\a)) + 
np(i\a) is its Doob decomposition as a martingale plus a predictable (in that 
case deterministic) increasing process. Moreover, if n > 1, and if X is an 
J-n-i measurable random variable, we compute 

E(Xe„(i)) = ^go(«)E«(Xe„(i)) 

a 

= Y,qo{a)Ea{X)p{i\a) = ^ E(XQ„_i(a))p(f |a). 

a a 

For the last equality we used the Q's characterization as Radon-Nikodym 
derivatives. This proves that E(e„(i)|J>i_i) = (5n-i(«)p(^|a) = "^nii)- 
Hence, setting 

n 

A(«) := E ^'"(^)' 

m=l 

an increasing predictable process, we find that X„(i) := Nn{i) — ^n(^) is an 
J>t-martingale under P, so each Nn{i) is again a sub-martingale with Doob 
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decomposition 

Nn{i) = Xn{i) + An{i) (19) 

under P. 

Finally, by some simple algebra we may rephrase the random recursion 
relation satisfied by the Q's as a stochastic difference equation 

Qn{a) - Qn-i{a) = Qn-i{a) J2 ^^(^n(i) - ^n-i(O)- (20) 

6.1.2 Derivation of the pointer state distribution evolution 

Equation (j20p admits eq. dlSp as a naive continuous time limit when 5, the 
scaling parameter, goes to 0~*". To put the validity of this formal approach on 
a firmer ground, one needs to prove the existence of a continuous time limit. 
This is a classical topic, but the presence of the scaling parameter S in various 
places prevents us from applying standard theorems straightforwardly. So we 
rely on a down-to-earth approach, which is rather technical. For this reason 
we relegated the argument to appendix [Bj This is where the interested reader 
should look for some background, precise definitions, etc. We give here a brief 
summary: 

- By an appropriate interpolation procedure, one defines a (5-dependent 
push-forward //a((5) of each Pq and fi{S) of P in Co(P^,R''^), the space of 
continuous functions from P"*" to P^ vanishing at 0. 

- We are not able to prove the convergence in law of the //q,(5) or of iJ,{6) 
when (5 — )• 0"*". 

- However, the finite dimensional distributions of the joint processes 
Nn{i) and Qn{ct) under each P, (where • stands either for an element of 
S or for nothing) admit, after appropriate time dependent centering and 
scaling, continuous time limits which are the joint finite dimensional distri- 
butions, under a probability measure /i, on Co(P~'', P^), for processes Wt{i), 
to be thought of as 

lim VS{Nt/s{i) -po{i)t/6), 

and Qtia), to be thought of as lim^_^o+ Qt/sict)- 

- The process Wt is the canonical coordinate process on C(P^, P''^), and 
its natural filtration Qt is to be thought of as the continuous time limit of 
the natural filtration for Nn, i.e. as the information collected by indirect 
measurements up to time t. 

- The identity = ^a1oi'^)f^a holds. The Radon-Nikodym derivative 
of /u(a) with respect to n on Qt is Mt{a)/Mt where 

M,{a) := eE«r(.Hm«-|E.Po«rWa)^^ :=^qoia)Mtia) 

a 
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For each a, Mj~ (a) is a ^(-martingale under /ia, and Mf is a ^(-martingale 
under //. 

- For each T > 0, under the measure M^^dfj, (which coincides with 
M^^ {a)d^a for every a), the process (W()(g[o,T] is a continuous time-homogeneous 
centered Gaussian process with covariance mm{t, s)[6ij po{i) — Po{i)po{j)) ■ 
Thus, by Girsanov's theorem, under each /Xq, Wt is a continuous time- 
homogeneous Gaussian process with independent increments, characterized 

by 

^^"{Wt{i)) = tpo{i)r{i\a) 
Cov^'"{Wt{i),Ws{j)) = mm{t,s){5ijpo{i)-po{i)po{j)). 

- There is an explicit formula for the Q^'s in terms of the VF^'s, namely: 

We are now in position to check that all the properties established in the 
discrete setting, as listed in section I6.1.H have a direct naive counterpart in 
the continuous time setting. 

The construction of the filtration Qt as the natural filtration for the 
canonical process was already explained. We have also already mentioned 
that there is an explicit formula for the Qt's. The counterpart of (fT8|) . the 
Laplace transform of the counting processes joint distributions is given for 
the canonical process in eq. (j29p . Appendix [Bl 

The counterpart of the counting process Doob decomposition under 
is Wt{i) = {Wt{i) — tpo{i)T {i\a)) + tpQ{i)T{i\a) under //q. 

To get the counterpart of the counting process Doob-Meyer decomposi- 
tion under P, i.e. the Doob-Meyer decomposition of Wt{i) under /i, we use 
Girsanov's theorem. As recalled above, for every T > (W^t(0)te[o,T] is a 
continuous martingale under M^^dfi. From 

dM(/M( = 5^go(a)^^5^r(i|a)dm« = J]Qt(a)^r(i|aW(i), 

a ^ i a i 

we infer that the increasing process 

At{i) := ds y^^Qs{a)po{i)T{i\a) 
Jo 

is the compensator of Wt{i), i.e. 

Xt{i) :=Wt{i)-Mi) 

is a Qt martingale under /x, with quadratic variation given by ()16p . It is easily 
seen that At{i), Xt{i) are the obvious continuous time limits of 
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It remains to write down the stochastic evolution equations for the Qt's. 
By Ito's formula for a ratio, we find 

Qt{a) V Mt{a) Mt ) \ Mt ) ' 
leading immediately to (fT5|) which we reproduce for convenience: 

dQt{a) = Qt{a) (r(f|a) - {Ti)t) dXt{i) 

i 

where (rj)^ := (3t(/3)r(i|/3). Note again that this equation is also the 
naive continuous time limit of the discrete equation (j20p . 

To summarize, one makes no mistakes if one works naively and forgets 
about the lengthy rigorous construction of the continuous time limit. This 
gives us confidence in what follows to proceed straightforwardly in the deriva- 
tion of continuous time equations in more complicated situations. 



6.1.3 Convergence of the continuous time evolution 

We now prove the convergence of Qt{oi) when t goes to infinity. Its almost 
sure convergence is a direct consequence of its martingale property. We need 
to prove that the final distribution is 

n ^ ifa^T 

^ " I '?o(a)/g-o(a) ifaGT ^^'^ 

and that the convergence is exponential with the characteristic time rxa. 

First we prove that the limit of the sector distribution Qt{,<^ '■= "^a&a Qt{(^) 

is 

Ooo(a) = 5a,T 
The equation of evolution for the sector distribution is 

dQM) = dQtia) = Qt{a) ^(r(i|a) - {r,)t)dXt{i) 

In the limit t — >• oo, we have Qoo(o)(r(i|a) — {Ti)cxi) = for all i, fi- 
almost surely. Then either Qoo(o) = or r(i|a) = X^j Ooo(/5)r(i|/3). Since 
r(i|a) 7^ r(i|/3) if a 7^ /3, the solution to the limit equation is Qooio.) = Sa,T- 
Second we show that = So(°Q. if q,' q, are in the same sector. As in 

Qt{a) qo(a) ' 

the discrete case, this relation implies eq. ()2ip . We compute: 

Qtjct) , Qtia') _ dQtja') _ dQtja) , dQt{a) .2 _ dQtja') dQtja) 
Qt{a') Qt{a) ~ Qt{a') Qt{a) ^ ^ Qt{a) ' Qt{a') Qt{a) 
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Since, '^'/".^^ = ^S^^"-) for a and a' in the same sector, we obtain (i^rr-r = 

' Qt{a') Qt(a) ' Qt(") 

if a, a' £ a. For all time t, 

Qt{a') qo{a') , 

= II O ~ Q 

Qt{a) qo{a) 

which achieves the proof for the limit pointer state distribution. 

Finally we prove the exponential convergence. The tools we use are the 
convergence of the pointer state distribution and the Ito calculus. From the 
distribution convergence, we have {Ti)t — T{i\T) for t large enough. The 
evolution equation ()15p for a T becomes dQt{a) ~ Qt{oi)Y^j^{T{i\a) — 
T{i\T))dXt{i). This equation is a well known stochastic exponential equa- 
tion. Thus, at large time t, with good approximation 

Qtia) ~ const. e-5E.PoW(r{i|a)-r{i|T))2+E,Xt(i)(r(i|a)-r(i|T)) 

for a T. Keeping only the leading term in the exponential we obtain the 
exponential decrease, Qtia) — exp(— t/rxa), with rxa given in eq. (fT7|l . 



6.1.4 DifTerent partial measurement methods 

The previous results can easily be extended to cases where different mea- 
surement methods are randomly used. Since proofs are similar to those of 
previous sections, here we only present a general outline of the approach. We 
limit ourselves to a time and realization independent partial measurement 
method distribution. In this case d„(oi,ii, • • • ,z„,o„+i) = Ilfcii c(ofc) with 

EoC(0) = l. 

To stay in the scope of the diffusive limit we assume that for any o, 
(il^-W) / 0. 

Following previous sections, we define linear interpolations W^^\o,i) of 
the counting processes which naively read 

V5{N,/s{o,i)-c{o)p"o{i)t/^) 

See Appendix ICl for precise definitions. As shown in this appendix, all finite 
dimensional distribution functions of W^^\o,i) under (a push-forward of) 
Pa (resp. P) have a finite limit as 5 — t- 0^ which coincide with those of 
continuous random processes, denoted Wt{o,i), under appropriate measures 
denoted fia (resp. fi). Under /Iq,, Wt{o,i) is a Gaussian process with 

E^'^{Wt{o,i) = tcio)p"Q{i}T^"\i\a) 
CoY>'''{Wtio,i),Ws{o',j)) = mm{t,s) (c(oK(i)<5(„,,),(,,,,) - c(oK(i)c(o'K'(j)) • 

with = and r(°)(f|a) = 2Im ("i^l^g^ 
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The measure is the sum /i = qQ{a)fia- The Radon-Nikodym deriva- 
tive of fi{a) with respect to n is Mt{a)/Mt where Mt = qQ{a)Mt{a) with 

Mt{a) = eS{°,«)eBr(°'{i|")i^t(o,0-|E(o,,)eB'=(°)PoWr(°){i|")' 
As in the section l6.1.2l we define 

Xt{o,i) = Wt{o,i) - fy^Qs{a)c{o)pl(i)T^°^(i\a)ds 
Jo „ 

The Xt{o, i) are martingales under From this definition we obtain straight- 
forwardly 

dQt{a) = Qt{a) ^ {T^°Hi\a) - {r^°Ki))t)dXt{o,i) 
(o,i)eE 

with (r(°)(i))t = Ear(°Hi|a)Qt(«) and 

dXt{o,i)dXt{o,j) = dt{c{o)po{i)5i^o,i),{o',j) - c{o)po{i)c{o')pQ {j)) 

The limit of Qt{ct) is the same but the sectors are now the sets of basis 
states such that T^°\i\a') = T^°\i\a) for all partial measurement methods 
and all partial measurement results. The convergence toward the limit dis- 
tribution is exponential 

Qi(a):^exp[-^ c{o)p°,{{){T^°\t\a) - r^'^\i\T)f 

{o,'i)eE 

The approximation hold if t is large enough. The convergence is exponential 
with a characteristic time 

— = E c{o)p°om("Hi\cy)-r^°\mf 

We find a convergence rate which is a mean convergence rate as in the discrete 
case. The same result is found by taking the scaling limit of the discrete case 
mean relative entropy. 



6.2 Density matrix evolution 

We are now interested in the density matrix evolution. 

As in section [5^ the density matrix at time n can be decomposed in the 
basis of pointer states: 

a,)3 
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The same decomposition applies to the time continuous density matrix we 
win define. The recurrence relation ([T|) translates for A„(a,/3) in 

^ , ^, ^„_i(a,/3)M(°")(i„|a)M(«")(i„|/3)' 
An[a,p) = — 



Where M^°\i\a) = {i\U^°^ {a)\^^°^) . For a = /3, this reproduces the pointer 
state distribution recurrence relation ([3]), as expected. 

We first limit ourselves to the case where only one partial measurement 
method is used and we omit the index o. The results will then be generalized 
to different partial measurement methods. We used a few hypotheses to get 
the continuous-time limit: 

- The first two assumptions are related to the development in of the 
conditional probabilities p{i\a). As stated before, the interaction Hamilto- 
nian must be rescaled Hj — ?■ '~^Hj and for any partial measurement result 
i, / 0. Then 

p{i\a) = po(^)(l + VSTs{i\a)) 

with po{i) = \{i\ip)\'^ . The assumption ^ leads to the diffusive limit. 

If this condition is not fulfilled for every i, then a jump-diffusion limit is 
found as shown in [23]. 

- A third assumption is needed to obtain a convergence of the evolution 
of the phases between different pointer states. The interaction Hamiltonian 
expectation must be zero : 

= 

Under these assumptions, we show in Appendix |D] that the time contin- 
uous evolution derived from the discrete time case is 

^l(a,p)t-ij:MiH-c{im*)Wt(i) 



with c{i\a) = ^^^fj^^'^'^ and 



l{a,p) := -i{E^ -Ep)-]- ^po«(|c(i|a)P + |c( 



2 , , c{i\a? - cmY'^^ 

i 



If we set a = /3 we recover the result on the pointer state distribution. 

A simple computation using Ito rules shows that this process is solution 
of a Belavkin diffusive equation: 

dpt = L{pt) - i Y^iCiPt - ptC\ - ptTr[{Q - Cl)pt])dXtii) (23) 
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with the Lindbladian 

L{p) = -i[Hs,p\ + Y^PomCipCj - \{Cja,p}) 

i 

and d := Ea c(«|a)|a>(a| = ^''^^f ^ ■ 

As shown in |20j . this equation corresponds to the time continuous hmit 
of repeated POVM processes ([1]) even if the non destruction assumption <^ 
is not fulfilled. 

In the next section we study the long time behavior of such evolution in 
the non destructive case. 



6.2.1 Long time convergence of the density matrix 

The pointer state distribution convergence indicates that, in the long time 
limit, the system is in a subspace of basis T. This information only tells 
us what is the limit of the elements At{a, f3) when a or f3 are not in the 
limit sector T. From the Cauchy-Schwarz theorem, limt^oo Qt{o)Qt{P) = 
implies limt^oo /3) = 0. For the elements At[a,P) with a,/3 G T, the 

limit t — )• oo is yet unknown. 

We decompose the operators Cj in a sum of two hermitian operators 

Ci = Ri + iSi 

with Ri = Ylia Re(c(i|a)) |a)(a| and Si = Ea |r(i|a)|a)(a|. 

As in the discrete time case, the density matrix evolution has to be 
modified by a unitary process in order to get convergence when t goes to 
infinity. Let Ut be the unitary diagonal operator defined via 

U^UUt = -i{Hs - Y,Po{i)[RiiSi - 2{Si)t) - \Rj])dt -iY^RidXtii) 

i i 

and let pt be the modified density matrix 

Pt = uJptUt 
As we show below it has an almost sure limit 

hm Pt = ^^Pt Po Pt (24) 

where Px := '^^^x |7)(7| the projector on the subspace corresponding to 
the sector T. Therefore, poo is equivalent to the density matrix we would 
have found if an initial Von Neumann measurement had been performed on 
the system. The unitary operator Ut only induces a rotation inside the limit 
subspace. 
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Recall that we only need to prove the convergence of the pt matrix el- 
ements corresponding to two pointer states in the same sector. From the 
Belavkin equation ()23p and using Ito rules, we find the evolution equation 
for Pt. 

dpt = '^Po{i){SiPtSi - ^{SiSi,pt})dt - i^{{Si,pt} - 2Tr[Sipt])dXt{i) 

i i 

Thus, the time evolution of matrix elements At{a,l3) of pt with /3 and a in 
the same sector is, 

dAt{a,P) = It(a,/3)^(r(i|a) - {r{i))t)dXtii) 

i 

Noticing that Qt{a)dAt{a, (3) = At{a, /3)dQt{a) and repeating the discussion 
of section I6.1.3| we get 

T . m f ifaand/3GT 
^-(«'^) = |(r^ else 

This proves the limit (|24p . 

6.2.2 Extension to difTerent partial measurement methods 

We can extend our results to cases where different partial measurement meth- 
ods are used. Once again we limit ourselves to time independent random 
protocols. The density matrix evolution is modified as follows: 

dpt = L{pt)dt+ ^ V(^o,i){pt)dXt{o,i) 
(o,i)eE 

with 

-Dio.M) = -^ict^Pt - PtCf^ - Pt Tr[cfpt - ptC\''^\) 
where ) = and 

L{pt) = -^[H,,p^]+ <o) ^{0^^ ptC[''^^ -\{cf ct\ Pt}) 

As before c(o) is the probability of using measurement method o. The limit 
density matrix can be analyzed as above: we obtain identical convergence 
statements once the density matrix has been rotated using an appropriate 
unitary Ut- 
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A Details for mutual singularity 



We prove that if o is recurrent but (o, i) is not then Yl^y Q{l)P° iAl) — 0- 
Observe that under we have the Markov property 

E«(l/„=i|J-„_i)=p«"(i|a) (25) 

It says that, under Pq,, 0„ = o the next measurement outcome is i G spec(o) 
with probabihty p°{i\a) independently of what happened before. 

Now assume that under Pq measurement method o is recurrent with 
probabihty 1. Take < Ti < r2 < • • • < Tfc • • • to be the times when the 
mcasiucment method is o. Wc show, using the strong Markov property, that 
It^ , , • ■ ■ are independent identically distributed random variables with 
distribution p°{-\a). This is quite natural: the functions d„ help choosing 
the measurement method, but they do not influence the measurement result. 

Indeed, note first that the above statement is trivial when there is only 
one measurement method, because then there is no need to invoke stopping 
times and the strong Markov property. In the general case, note the slight 
mismatch with usual notations: {T^ < n\ is in fact Fn-i measurable, so it 
is natural to write Tt^-i for the algebra associated to the stopping time Tj.. 
Then write 



= l7Ti=n • • • '^iT^_^=ik-iP"{'^k\a) 

One can go on to condition with respect to Ft^_^-1i ■ ■ ■ until one finds the 
plain expectation 

Ea(l/Ti=n • • • 1/Tfc=ifc) =P''{h\oi) ■ ■■p°{ik\a). 

As a consequence, for any a such that o is recurrent under P^: 

- either p°{i\a) > and with Pa-probability 1 the outcome i appears 
infinitely many times in the sequence It^, 1x2^ " ^ (o, i) is recurrent with 
probability 1, 

- or p°{i\a) = and i never appears in the sequence It^^It^, ■ ■ ■ , i-e. 
(o, i) never appears. 

Now assume that o is recurrent under all Pq's. The above implies imme- 
diately that the probability under P that {o,i) is non-recurrent (this event 

is denoted by i(o^j)) is given by J2-y,po(i\^)=oQo{l)- 

If p"{i\f3) = then Q{P) = Q{p)l^ because by the recursion relation 

Qn{P) = whenever (o, i) has shown up before time n. So 

E(g(,)|i,,,)-^(^(^))- -(^) 



E(1a,„,J E^,po(^|^)=o5o(7) 
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which imphes that 

E( Yl Q(7)|i(o,.)) = 1 

7,p°(i|7)=0 

Hence, conditional on ^{o,j)i the Q{a)^s for which p"{i\a) > have to vanish. 
Equivalently, Q{'y)p°{i\'y) = for each 7 and Q(7)p°(^|7) = 0, which was 
to be proved. 

B Proof of existence of a continuous time limit 

We first put the notion of continuous time hmit in context. 

Let V be the vector space Co{R~^, K^) of continuous functions /. from R+ 
to such that /o = G R^. For each 6 > 0, and each a; G we define a 
function wl:^\i) on by Unear interpolation of wl:^\i) := \/6{N^/g(i) — 
Po{i)t/5) if t/6 is an integer. Explicitly, for t G [6n,6{n + 1)] 

Wt^^\i) = V6{{n + l- t/5)Nn{i) + {t/5 - n)iV„+i(i) - po(i)i/5) • 

For every uj ^ Q. the function W^^\i) is continuous for t G R"*". So we 
have a map l^^"^) : — )- 1/. But, as already pointed out before, there is no 
hope that, for a fixed a; G fi, Wi^\i) has a limit when 5 — )• 0"*". The only 
clear fact is that for a fixed i, the central limit theorem ensures that the 
distribution of W^^\i) under each Pq, has a Gaussian limit when 5 — )• 0"*". 
Note that this observation fixes the scaling \/^ as the only one possible. 

But if we are interested in convergence as a process, a deeper approach 
is needed. If we endow V with the topology of uniform convergence on 
compact sets TiV) and with the corresponding Borel cj-algebra B{V), we 
can show that the map W'^^'^ is measurable from to {y^B{V)). This 

is not difficult, because by a classical result, B{y) is the smallest cr-algebra 
on V containing the family of sets 

Bt,^,a :={/eCo(R+,R^),a</t(i)} 

indexed by t g]0,+oo[, i G / and a G R. It is plain that the inverse image 
of 5*'*'" under H^(^) is in whenever n > t/5. As the filtration on O is 
exactly the one making the Nn J^n-measurable, the appropriate filtration 
on C(R+, R^) should be the natural one, the smallest making the canonical 
process adaptec^. We denote it by Qt- 

Then any probability measure on (^2,7^) induces via I^('') a probability 
measure on (y,B(y)). Note that, via (fT^ . the measures we defined previ- 
ously on [Vt) depend on 6, and to be explicit we write rW, , etc. to stress 
this fact. Let jiai^) be the image measure of Rq ^ pushed forward by W^^^ on 

®Of course, as long as 5 > 0, Wt looks a bit forward in the future, as it involves Nn+i 
for t £ [5n, 5{n + 1)], but for continuous process this does not matter. 
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{V,B(y)). As {V, T(y)) is a so-called Polish space, there is a nice notion of 
convergence for measures on it, called "weak convergence" of measures, and 
we could ask if the ^q,((5)'s converge weakly to some probability measure on 
{V,B{V)) (and then, so would fi{6) = X^q, '?o(ct)/"a('^))- Note that despite its 
name, weak convergence is strong enough to ensure the convergence of the 
expectations of rather general functionals, so we could hope to control the 
Q's continuous time limit as well, because they are nice functionals of the 
counting process. 

It is usually in this context that continuous time limits have a meaning. 
In this setting, there are a number of theorems, called functional central limit 
theorems, or Donsker invariance principles, that express the continuous time 
limit of random walks (with independent increments) in terms of Brownian 
motions in fine details. Alas, though under each P^, Nn{i) is a random 
walk with independent increments, the theorems we are aware do not apply 
immediately. The problem is that the 5 dependence is not only in W, but 
also in the P^'s. This problems would show up even more dramatically to 
deal with the Q's convergence and there relationships with W. While we 
think these are purely technical details in our case, we shall not try to deal 
with them: we shall instead rely on a weaker and slightly less natural notion 
of continuous time limit that will suffice for our purposes. To say things in 
more mathematical terms: we shall content ourselves with proving that the 
joint finite dimensional distributions of Ws and Q's converge to joint finite 
dimensional distributions of continuous processes we can identify explicitly, 
but we do not embark on the more technical task of proving tightness. 

We now turn to our explicit approach of the continuous time limit. 
From the characteristic function ()18p we obtain easily that 

lim (e'^'-^^^^'^'^^<^'^-<-^'^^ 

\ a 

elEti(<!-ti-i)(E,poWAiW'-(E.po«A,«)') (26) 

for /c > 1, arbitrary non-decreasing sequences = to ^ ^ " " " ^ ^fc of 
length k and arbitrary (complex) A/(i)'s. 

By a standard theorem on characteristic functions, we have thus proved 
that the wl^\i) finite marginals (under each P^^ and under P^*^^) have a 
limit for 5 — )• 0"*". This is much weaker than what weak convergence of 
measures would ensure. It is enough to ensure that the limit marginals satisfy 
the Kolmogorov consistency criterion, but it does not guarantee that it is 
possible to concentrate the corresponding process on Co(R^,R''). However 
in the case at hand, we can bypass this problem because of the simple form 
of the result, which is Gaussian for each a. 
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Let V be the Wiener measure of a standard Brownian motion on Co([R^, K^). 
The hnear map from to itself defined by 

3 

induces a map from Co(IR^, K^) to itself. Let be the image measure of v 
under this map. It is easily seen that under this law, the canonical process 
on Co([R+, K^) satisfies 



for A; > 1, arbitrary non-decreasing sequences 
length k and arbitrary (complex) Xi(iys. 

As a consequence, Ito's formula holds with 



^0 ^ ^1 ^ ■ ■ ■ ^ ^fc of 



dWtii)dWtij) = dt {5i,jPo{i)-poii)poij)). 

Our aim is to use Girsanov's theorem to deform the measure and go 
from the right-hand side in (|27p to the right-hand in (12611 . If the process Ut 
with values in is adapted and satisfies some further technical integrability 
conditions, 

Mf := ei'oS»^=W'^^=«-|/o(E,poWC/.W'-(E,poWC/.W)')'i^ 

is a martingale. Moreover, by Ito's formula, dMj^ = M^^ '^■Us{i)dWs{i). 
Then, by Girsanov's theorem, for any T > 0, under the measure dfij, : = 
M^dfi"^ on Co([0, T], R^), the process Wt{i)-po{i) Jo Po{j)Us{j))ds, 
t G [0, T], has the same law as Wt{i), t G [0,T], under dfi^. 
Note that for t < T, Mj^ is a Radon-Nikodym derivative. 



Mr-- 



dfij, 
~d{^ 



( d^j, 



d^i^ 



In general this construction cannot work for infinite T, because /Xji and 
become singular. However, T plays only a dummy role : for T' < T, 

dfij^ and d/i^, coincide on Qt'- So it is only a slight abuse, which lightens 

notations a bit, to write d/i^ for d/x^ and 



u 



dfj, 



dflO 



dfi'- 



dfi^ 



For the special choice Ut{i) 
compute 



r(i|a), using Yl^iPo{i)T{i\a) = 0, we 



Mt{a) := e 



E,r(i|a)w^t(i)-|E.po(i)r(i|a)2 
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which is certainly a martingale, and we obtain for every T a measure dfj,a 
on Co([0,T],lR^) such that 

Er=i(t!-ii-i)E.A,(i)poWr(i|a)giEti(*i-ti-i)(E,poWAiW'-(E»po«AKi))') (28) 



e ' 



for k > 1, arbitrary non-decreasing sequences = to ^ ^ " " " ^ ^fc ^ ^ of 
length k and arbitrary (complex) A;(i)'s. 

Finally, setting Mj := '^^qQ{a)Mt{a) (trivially a martingale again) and 
diJ, := Mtdjj^ we obtain for every T a measure d/Li on Co([0,T], R^) such that 



EM (^eEtiE,A,W(M/t,»-m,_i(i))^ 



Etife-ti-i)E,A,»PoWr(i|a) 



elEti(<!-*i-i)(E,poWAiW'-(E,po«A;(i))2) (29) 

for fc > 1, arbitrary non-decreasing sequences = to < ti < ■ ■ ■ < < T oi 
length k and arbitrary (complex) A;(i)'s. So we have found the continuous 
time limit of the counting process as the canonical process on Co([0,T], R^) 
with the measure fi. 

It remains to deal with the continuous time limit of the Q„'s. 



Mtja) 



We note that by the chain rule ^ 



dfi{a) 
dfj, 



Qt ) is the 



Radon-Nikodym derivative of /i(a) with respect to /x. So with some mem- 
ory of what happened in the discrete case, it is natural to define Qt{a) : = 

/ N Mt (a) 

Using the explicit formula for the QnS in terms of the counting processes, 
we define Q^^'' by an interpolation procedure: 

Ql \a) := qo{a) jr-, , 

E/3 QoiP) U^ Pm"^^ \i)/V-S+Poii)t/6 

so that if t/6 = n, an integer, Q[^\a) = Qn{a). Note that in this formula 
the p(i|/3)'s depend implicitly on 6 via (jl4p . 

One can prove that the joint finite dimensional distributions of the pro- 
cesses {W^^\q[^^) under P^*^^ have limits when 5 — )• 0^ and that these limit 
are nothing but the joint finite dimensional distributions of the processes 
{Wt,Qt) under /i. In this precise sense, the QnS continuous time limit is 
deciphered. 

This result is really no big surprise, but to prove it we have to rely on 
an ad hoc trick and an explicit elementary but tedious computation. The 
details are neither illuminating nor elegant so we omit them. 
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When 6 is small enough, all non-empty sets in J^n have strictly positive 
measure, so that if qo^a) > the same is true for Qnice) for all n's. As 
furthermore Qn{ct) = 1 for all n, all the information on the QnS (joint) 
laws is embodied in the joint laws of ratios of QnS. 

As these ratios have a simple product structure in terms of the counting 
processes, the explicit computation of 



gEf^l Mi){Nni (i)-^n;_i (i)) -Q -Q 

(=1 a,/3 



for A: > 1, arbitrary non-decreasing sequences of integers = uq < ni < 
• • • < rifc of length k, and arbitrary (complex) A/(i)'s and r/i(a,/3)'s is in 
some sense a special case of (jlSp . 

The same remark applies to computations involving ratios of Qts. 

This allows to compute explicitly that 

/ k / (S) \ ^i('^^l^)^ 



5^0+ 



1=1 a, 13 



£M I^^Et, E. hUmt, ii)-w,, , (.)) ]J J| j^ ^*' ^ "^"'"^^ j (30) 

for /c > 1, arbitrary non-decreasing sequences = to ^ ^ • • • ^ ^fc of 
length k, and arbitrary (complex) Ai(z)'s and r/i(a,/3)'s. 

As such a mixture of Laplace and Mellin transforms characterizes the 
distributions completely, this concludes the existence of a natural continuous 
time limit. 



C Details on the continuous time limit with differ- 
ent partial measurement methods 

We use the linear interpolation of appendix[B]on W^^'' (o, i) := \^{Ni/sio, i) — 
c{o)pQ{i)t / 6) if t/6 is an integer. Explicitly for t £ [6n,6{n + 1)], 

wi^\o, i) = V6{in + l- t/d)Nnio, i) + {t/5 - n)Nn+i{o, i) - c{o)pl{^t/5) . 

We remind that E = \Jg^(^{o} 0spec{o) is the set of all possible measure- 
ment methods and outcomes. We expect the limit time continuous process 
to live on the vector space of continuous function from R"'" to . 

Compare to previous sections, the main changes are in the correlation 
functions calculations. Thanks to the measurement method distribution time 
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and realization independency, we find: 

k / \ '''""'-1 

E^o(«)n Yl e''^^''^cio)p%z\a)\ (31) 

a 1=1 \(o,i)eE / 

for /c > 1, arbitrary non-decreasing sequences of integers = uq < ni < 
• • • < rife of length k, and arbitrary (complex) A;(o, z)'s. 
As a consequence, in the limit (5 — >■ 0+, 

lim E(^) (Qj:f=i^(o4)eEMo,i){Wt^P{o,i)-wip_^{o,i)) 

^^gQ(«)eEti(*i-*i-i)E(o.i)eBAKo,i)c(oK(i)rW(i|a)j ^ 
elEf=i(t(-t!-i){E(o,i)eB'={°K«A((o,i)2-(E(o,i)eBc{oK«AKo,i))2) ^32) 
for > 1, arbitrary non-decreasing sequences = to < ti < ■ ■ ■ < tk of length 



k and arbitrary (complex) A/(o, i)'s and with T^"\i\a) = 2Im 



Then each Tyj'^''(o, i) under /j,^^^ converges toward a process Wt{o,i) under 
fi. 

The demonstration is then as in previous section except for notational 
differences which keep track of the o-dependency of Wt{o,i). As in the 
previous section, the measure fi^ is defined as the push-forward measure of 
u on under the linear map 

y{o,i) := c{o)pl{{){x{o,i) - c{o)p^{i) ^ c{o')p^ (j) x{o' , j) 

{o',j)eE 

The martingale Mt is defined by 

Mt = Y^ go(c)e^(''.i)6^5r('''{i|")w^t(o,0-| E(o,i)eBc(oK(i)r('')(i|a)2 

a 

The measure fi is defined via Girsanov's transformation: E^(-) = E'*°(Mf). 

D Derivation of the density matrix evolution 

Let us derive the density matrix continuous time limit. Recall that at time 
n its elements are functions of the counting processes 

Ao{a,f^)U^{M{^\a)Mm*)^-(^) 
'^^ E.^o(7)n.p(^l7)^"« 
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from this expression we define time continuous processes 



,(5). a, _ Ao{a,^)lli{M{i\a)M{ii 



E,^o(7)n.p(^i7)^* 

equal to A„(a,/3) if t/6 = n is an integer. The M(i|a)'s depend explicitly 
on 5 via 

M{i\a) = (f|e~^'^^^"'+^''+75^-)|^) 

We rewrite the products over the partial measurement results as exponentials 
of sums 

4 )(a,/3) = °^ '^^ 

^^gQ(^)eEiln(p(i|7)/P0»)(w/''»/v^+P0»|) 

A detailed analysis of the limit (5 — >• 0"*" would require to perform the same 
study as in section [Bj However, at this stage we are confident enough to state 
that we can safely shortcut a few steps and use directly that W^^'^ (i) converge. 

Using {^\Hi\^) = 0, and the identity J2iPo{i)^^^^ = ZiPo{i)\c{i\a)\\ 
we obtain 

hm Vln(M(i|a)M(f|/3)Vpo«)(W^f^«/^^ + Po(i)|) 

i 

= l{a, p)t - i EMi\a) - c{i\mWt{i) 

where the limit as to be understood as the limit of any finite dimensional 
correlation functions. Therefore 

lim^f^(a,/3) =At{a,^) 

with At{a, P) defined in (|22p . It is then a simple matter, using Ito rules 
for Wt{i), to derive the Belavkin equation (j23p for the density matrix pt = 
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